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By liax A. Haaelet 

SUMMARY 


In this report, graphs of constants are given which 
are to he need together with a knowledge of maximum total- 
head Iobb , static— pressure decrement, and the Integral of 
total— head loss across the wake of an airfoil for the rapid 
determination of section profile— drag coefficient. The 
constants './ere computed under the assumption that the total- 
head pressv.re loss In the wake of the airfoil has a cosine- 
squared distribution and that no variation in static pres- 
sure exists across the wake at any given position. The 
range of pressure losses, for which the results are given. 

Is sufficiently large for usual wind-tunnel and free— flight 
data, r.nd conpress lb lllty effects are considered for Hacli 
numbers up to and Including M ■ 1, where M is the Ilach 
number of the free stream. 

Included In the report are results of computations 
that were carried out, for free— stream Ilach numbers between 
0.5 and 1.0 and for certain assumed types of total— pressure 
distribution in the wake, to compare theoretical drag 
coefficients as determined by various equations based upon 
the momentum method. Among the assumed wake shapes are 
forms similar to those encountered at supercritical Bpeeds. 
Results obtained by point— by— point methods of Integration 
are compared with those computed by means of the constants 
mentioned above. Yor the oases examined, the numerical 
agreement becomes less satisfactory as the maximum total- 
head decrement increases but for such total— head losses as 
are usually encountered In practloe the agreement is quite 
sufficient. It is therefore concluded that, within the lim- 
its of the validity of the basic assumptions underlying the 
momentum method and for total— head decrements of normal 
magnitude, the use of the more rapid technique Is Justified 
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for vide variations in wake shapes for the Kach number 
range considered. 


IHTSODUOT IOIT 


One accepted method of determining the profile drag 
of an airfoil is developed from consideration of the 
momentum defect in the' wake of the body. As given by 
Jones (reference 1), for velocities at vhioh air may be 
assumed an incompressible medium, the actual computation 
involves the evaluation of an integral with integrand a 
function of total— head and stat ic— pressure losses across 
the wake. In order to achieve this evaluation in terms of 
measured quantities, it is customary to presuppose that 
there is no mixing in the vake in stream tubes between the 
plans of measurement and a second plane vhioh is far enough 
dow.istroam that the static pressure there may be assumed 
to have returned to its original free— stream value. 

!fhen greater velocities are involved, it has been 
found possible to modify the Jones equation for compressi- 
bility effects so that the desired section drag may once 
more be found in terms of the same preBBure measurements. 

As a counterpart to the previous assumption concerning 
mixing, it is again supposed that streamlines may be drawn 
in the wr.?"9 , between the two planes, so that total head 
is constant along oach stream tube. An added condition 
is necessary, however, in order that density may be eval- 
uated and, to this end, the total energy per unit hibbs of 
air is assumed constant in any section across the field 
of flow. 

The acceptance of the momentum method as a valid way 
of finding the drag of an airfoil eection is dependent 
more on the experimental evidence at hand than on the un— 
assailabllity of its underlying principles. In reference 
2, G. X. fp.ylor has shown theoretically that, for oertain 
types of pressure distributions, the indicated drag may 
he as much as 10 peroent in error depending on whether the 
wake downstream mixes or flows in a streamline manner. In 
references 3 and 4, however, the method appears to have 
been substantiated experimentally by measurements taken at 
different distances behind the airfoil and by comparisons 
with values measured on a balance. The conclusions pre- 
sented in reference 4 indicate that the values of drag 
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coefficient obtained by means of the momentum method are 
reliable up to speeds near the velocity of sound and that 
the presence of shook waves of limited extent does not in- 
validate tlie result. These conclusions must not be con- 
sidered as a complete negation of Taylor's results, since 
his investigation of the pressure distributions encountered 
in practice, as given in reference 1, shoved that the error 
due to mixing vas in the worst case much smaller than was 
theoretically possible. 

The expressions for drag coefficients, for either the 
compressible or incompressible case, appear as integrals 
which are to be evaluated in the measurement plans. Ac- 
tual values of head loss and static— pressure variation can 
be determined by means of a pitot— traverse method and, as 
a oon.Becv.ence, the given integrals can be evaluated by 
numerical integration. In order to avoid the lengthy com- 
putation Involved, an alternate approach has been developed 
by Sllversteln and Katsoff (reference 6) in whloh it is 
assumed that drag is proportional to the Integral of the 
total— head loss acroae the wake. The integral of head loss 
is obtalnod from an averaging rake, an integrating manometer, 
or by an integration of measurements made by a traversing 
tube. The proportionality factor, in turn, is evaluated 
by assuming a definite type of wake shape and, in this re- 
port, is tabulated as a function of maximum total— head loss, 
static— pressure decrement, and Ilach number of the free 
stream. The results given herein deviate somewhat from 
those given in reference 5 for reasons that will be dis- 
cussed later . 

In reference 6 a comparison was made between values 
of drag coefficient computed by means of point-by— point 
integration and those determined by the Integrating method 
and proportionality constants of reference 5. It was 
found that the latter procedure gave results in excellent 
agreement with those obtained by the former, at least up 
to values of the Mach number of the free stream in the 
neighborhood of 0.6. A vide range of shapes was consid- 
ered and, remarkably enough, the agreement remained uni- 
formly good for the Mach numbers considered. 

Since in high-speed wind tunnels it is necessary to 
evaluate drag at rlach numbers above those considered in 
reference 6, and since the recalculation of the nropor— 
tionallty constants resulted in a change in their previ- 
ously determined values, the present paper has undertaken 
the comparison of theoretical drag coefficients obtained 
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by differ ont methods for the high-speed range. She ap- 
pearance of shbck waves on an airfoil may-be expected at 
some value of free— stream Haoh number less than unity and, 
as a consequence, the distribution of total— head loss 
differs considerably from that usually encountered. Some 
of the \ro.l:o shapes considered are therefore chosen to 
correspond approximately to this situation. 


0UTLIF3 07 H3TH0D 


Complete details concerning the development of formu- 
las for the drag coefficient may be found in the refer eacee. 

It is considered sufficient for the present purposes merely 
to list the main results. The symbols used are defined in 
the appondix. 

In figure 1, the airfoil section is indicated together 
with the three planes of eBpeclal importance to the theory. 

Plane 0 is far enough upstream that free— stream conditions 
may be assumed to exist uniformly across this section. Plane 

±y 

Prom momentum considerations the drag d, per unit _ 
length of the airfoil, is 


3 is at a .'-rent enough distance downstream that the value 
<?f a ta t-l.c-.-'jr.Q.asur e has returned to its original .value ! p 0 
Plane 1 is the actual plane of measurement. ** 


d = jf p 8 V a (V 0 - * a )dy 8 (1) 

+■ d* f A X (So- v * «/ t 


cfer Vt v* Ml TT^Jh d f /*, K - 

where tho integration across the wake is in plane 
proof of this relation, valid for Bubsonic compressible 
flow with possible limited shock waves, is given in ref- 
erence 4. 

The profile— drag coefficient for the airfoil section 
is an immediate consequence of equation (l), so that 
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If the fluid is assumed to he lnoompr eanihift_ awd continu- 
ous, and If no mixing .la supposed to take _ place In the 
stream tubes between planes 1 and 2, equation (2) can be 
written as ' ' W) 

( y \ h - v ; 


2 n / .P - P H — H I / H-H i 

c a n - / /l 1 —2. S_ 1 I 1 - / 1 2_ 1 d 7l (j 

° J v * H 0 — p 0 H 0 — p 0 ** H 0 P 0 p_v. 

*/• 

From this lorn of the equation, the drag coefficient is 
directly calculable once the distributions of sAat lc— 
nresaure and total— head loss across the wake a re known. . 

If It 1 b assumed, ae Silver stein and Katzoff did in 
reference 5,' that the Integral of total— head loss across 1 
the wa’:s is proportional to the section— drag coefficient, 
a proportionality constant 7^ may be introduced so that 

•* ■ s“ 7 1 5 (r^) « 

c J W H o“*P 0 c ^o-Po^y 


H-H 

1- d 7l (3) 

■“o — Po P- V. 


where the subscript av . denotes average value across the 
wake. In roforences’ 5 and 6, is tabulated as computed 

Pi-P ( 

from equations (3) and (4) under the supposition that 

H 0 — Hi H 0 -p ( 

is constant and that — — has a cosine— squared die— 

H o-Po 

trlbution; that is, 


H o~ H x 

H o“Po 


H 0- H x 

H o->o 


where |y| = — and is measured from the oenter line of 
2 

the wake. In figures 2(a) and 2(b) extended values of P, 


are given in graphical form as functions of 
/ H„- H.\ 


H o-Vuax 


P — P 
P( 
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Wlian compreBS ib ility effects are to he considered, 
factors jT c or J 0 /l ^ . are used, so that 


„ Ls 


,r izii a,, - £» f> I ( 5) ^ 
■K *6-*o- *i L- 1 J 




The evaluation of 
oeeds as follows: 


J 


F 0 £s sim^4r. to that Of P'j and’ pro— 

Equation (2) is rewritten In euoh a form 
that the integrand may he calculated directly from the ■ 
free— stream Ilaoh number, the static— pressure loss, and the 
total— head Iosb in the wake. .If oonstant values of the 
first two parameters and a cosine— squared distribution of 


the latter are assumed, it is possible to find 


and 


subsequently to solve equation (3) for T Q « In figures 
3(a) to 3(e) will he found T c /T j plotted as a function 


of K, 


(-“rrO ■ 

p(j/ 


and SJLZIA, 


E o ~ Po 


In this report, two forms of equation (2) are used 
in the point— by— point integration for the compressible 
case. The first expression. 



is derived in reference 4 under the assumptions that total 
energy per unit mass of air is oonstant in any section 
across the field of flow and that between planes (l) and 
(2) Bernov.lli's equation for oompresslble flow holds. 

This form has a great advantage in that auxiliary tables, 
which are Included in the reference, may be prepared so 
that the computations Involved are expedited. 

Ab dorived in reference 6, the alternate form of the 
expression for drag coefficient is 
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The expression 1 + ti , which appears in equation (7), is 
the so-called compressibility faotor and is a function of 
Mach number , 


1+ ti 


1 + a! + 

4 



K a 

1600 


M' 


80000 


Figure 4 shows this function, for values of M up to 1, 
in a form that can be used in numerical integration. The 
subscripts on n , as given In equation (7), Indicate 
the plane in which M is evaluated. 

One outstanding simplification possible in the inte- 
gration of aquation (7) was made by Bay H. Ifrlght of the 
Langley Memorial Aeronautical Laboratory In some unpub- 
lished results. If the second radical, in the Integrand 
of this equation, is written in t he form 


-°( 


x -s»> 
*0 J 


it may bn shown that 


is primarily a function of Mach 


number and secondarily dependent on 


1 - iS. 


This follows 


from the definition of 0, 


for 
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a a 




o-s) 

• \‘»V7 


A graphical representation of 0 Is given in figure 5. 

fhen proceeding to the actual integration, it is suf— 

fioient to haow the variation of -2 1 . the Haoh num- 

H o ” Po _ _ _ 

her of tho free stream, and the value of S a. —S E. In the 

H — u 
o *o 

computation it will he found necessary to use the indi- 
cated "ir os cure ratio)/ — ^2-, Tigure 6 gives values of 

A p 0 

this parameter as a function of free— stream Mach number, 

K. 




As mentioned previously, the values of P 0 /Pj. In 

figures 3(a) to 3(e) differ in most caaeB from the corre- 
sponding values given in reference 5. Since the original 
assumptions are the same, there seems little douht that 
this discrepancy is produced hy differences in the manner 
in which the Integration of the haslo equation is per- 
formed. This report assumes a distribution for total— head 

H — H 

decrement ; while reference 5 assumes a distribution 

E o“Po V , 0 8 -V a 8 

for velocity— squared decrement v . .... ■ „ t although the pre— 

V s 

*o 

else nature of this latter assumption is not stated ex- 
plicitly. An analysis of the two sets of results giveB 
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credence to the conclus ion that Sllverstein and Eatsoff 
have set 

iziS* . lalzlal 

Ho-Po V 0 » 

whloh is In agreement with the theory of incompressible 
fluids. Such a relationship does not Introduce large 
errors unless the Mach number Is of considerable magni- 
tude. but in this latter oase the relation Is untenable 
and the values of 7 0 /*i In this report would appear to 

be the true values. (These remarks apply -equally well to 
reference 6, since the graphs of P 0 /Pi given there are 

In agreement with the tabular data In reference 5. 


COUPTITAT IONS 


By means of the graphical data that have been pre- 
sented, and the tables Included in reference 4, It le now 
possible to determine values of drag coefficient from 
equations (5), (6), and (7) and known wake distributions. 
Such data can, in turn, be used for a comparison of the 
results obtained from point— by— point methods of Integra- 
tion and the Silverstein and Eatsoff Integrating method. 

In the five oases listed below, an arbitrary distri— 

but ion for ~ - ° -1 is assumed, together with fixed val— 

H o — Po 

p J — p 0 

ues for — . The drag coefficient Is then computed 

Ho “ Po 

for different values of free— stream Mach number. In 
ohoosing the wake forms, all of which are shown in fig- 
ures 7(a) and 7(b), an attempt was made in cases I, III, 
and 17 to present distributions that are highly distorted 
from usual distributions j while in cases II and 7 shapes 
were ohosen that resembled somewhat the types encountered 
when shock waves are present on the surfaoe of the airfoil. 

It will be noted that all the oases are concerned 
with symmetrical configurations. Typical asymmetrical 
shapes could be constructed by Joining two halves of cosino 
squared distributions with equal heights and different base 
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widths. Trom the theory underlying the derivation of con- 
stants , and J c /l| , It follows that for suoh typical 
asymmetrical shapes the values of drag coefficient given hy 
equations (4) and (6) will always he exact. Other forms of 
asymmetry may he thought of as deviations from these hybrid 
cosine— squared curves and the study of the variation in 
drag coefficient produced hy these deviations can he related 
directly to problems suoh as are considered here. 


Case I 

As an extreme example, the distribution of head loss 
is made rectangular. ThuB , as shown In figure 7(a), let 

H — H 

— = 0.10 for | y | ^ a and equal to zero for all 

H o -p o Pi-Po 

other values of y; moreover, let - — ■■ « 0.1. It fol— 

H o-*o 

lows Immediately that 

X r 4 y a 0.20 

° ' w a o ~ Po 

In table I, results of the computations are given. 

3cuations (6) and (7) should give the same numerical 
result sine? they are equivalent. The deviations that ap- 
pear are probably attributable for the most part to small 
errors arising In reading the graphs and, to a smaller ex- 
tent, to the -Interpolation necessary In the tables. The 
differences in results from equations (6) and (6) are prob- 
ably due to the variation between the actual wake form and 
the hypothetical shape used to compute although 

numerical inaoouraoles must again be considered. It is 
observed, however, that the percentage of error Is small. 


Oase II 

The distribution of head loss for this case 1 b shown 
In figure 7(a). Under the assumption that 

p , — p ft 

— * *• = 0.1, the drag coefficient has been computed by 

®o Po 

means of equations (5) and (7) for different Mach numbers. 
The results of these calculations are given in table II 
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and show that In each Instance the error In the data deter- 
-mlned by equation .( 5 ) .. i s„-small-. 


Case III 


Assume now that 


H e ~ Hi 

H o - Po 


0.30 for 



as 


shown in figure 7(a). In table III, drag coefficients are 
given corresponding to this distribution of head loss when 


Pi I Po 
Hq “ Po 


0.1, 0.0, -0.2, -0.4, respectively, and for vari 


oub Mach numbers. The calculations are based on equations 
(5) and (6). The disparities between corresponding answers 
are, in general, larger for these examples than in case I. 

On the other hand, the percentage of error seems relatively 
small in comparison with the deviation between the true and 
assumed wake shape. 


Case IV 


In this case, drag coefficient is computed from eaua- 
tions (5) and (6) when M * 0.4, 0.6, 0.8, and 1.0 under 


the assumption that 


Hq - H, 

Ho ~ Po 


0.10 for 



as shown 


in figure 7(b). Table IV lists the results for 

Ei s 0.0, -0.1, -0.2, and -0.4, and in each case the 

H o “ Po 

percentage of error Is small. 


Case V 


and for 


7or the distribution of 
Pi " Po 


H -w H 

— 0 — — i. shown in figure 7(b) 


H 0 - 


H o - Po 

0, -0.2, and -0.4, drag coefficient is 


computed at free stream Mach numbers equal to 0.4, 0.6, 0.8, 
and 1.0. In this case the maximum total-head loss is 0.7 
and, as can be seen in table V, the percentage of error is 

Pi - Po 
H 0 - Po’ 


quite large, especially for higher values of 
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CONCLUSIONS 


The computations carried out in this report show that 
the Integrating method of Sllversteln and Katsoff gives re- 
sults In good agreement with the values of drag coefficient 
computed by the much more laborious process of point-by- 
polnt Integration so long aB the maximum total-head loss is 
not too large. The limitation on the total-head loss Is 
necessarily dependent on the accuracy desired, but for such 
values as are usually encountered in practice the cases 
considered do not show exorbitant deviations in the results. 

The agreement in results holds for considerable varia- 
tion of the wake form and. as far as theoretical resultB are 
concerned, extends up to a free-stream Mach number equal to 
1. Since previous reports have Indicated that the momentum 
method Is not necessarily nullified by the existence of 
shock waves of limited extent, it follows that the Integrat- 
ing method need not be Invalid beyond the critical speed of 
the airfoil. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPUKD IX 
SYMBOLS 


7 velocity 

p static pressure 

q dynamic pressure 

H total pressure 

p density 

c chord length of airfoil 

7 distance measured vertically from wake center 

w wal:e width 

K Inch number, ratio of the stream velocity to the 
■ local velocity of sound 

d section profile drag 

c^ section profile— drag coefficient 

7 C pr o-jor t ionality constant for compress ibis flow 

proportionality constant for incompressible flow 

1+n conroressibility correction factor, q ■ a ~ = 

H K 2 l+n 

y rrtlo of specific heat at constant pressure to spe- 
cific heat at constant volume 

Subscripts o, i, 8, refer to conditions existing in 
three different planes. Subscript o denotes free— stream 
conditions, x denotes the plane of measurement, and a 
denotes the plane aft of the airfoil in which the static 
pressure has returned to Its free— stream value. 
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TABLE I.- ca AS A FU1TCTI0I OP MACH HUMBER 
FOR WAKE POEM BT CASE I 


Mach 

number 

c a 

( equation 

(7)> 

c a 

(equation 

(6)) 

c a 

(equation 

(5)) 

Percentate of error in 
equation (5) as compare! 
with equation (6) 

0.5 

0.167S 

O.167S 

O.16S7 

0.5 

.6 

.l6l4 

.1616 

.1624 

' .5 

•7 

.1550 

.1551 

.1556 

•3 

.3 

.1474 

.1476 

.1403 

• 5 

• 9 

,i4oo 

.1403 

.1409 

.4 

1.0 

.1334 

.1331 

•1333 

.2 


TABLE II.- c d AS A FUITCTIOH OF MACH HUMBER 


FOR WAKE FORK Ii'T CASE II 


Mach 

c d 

c a 

Percentage of error in 

number 

(equation 

"(7)) 

(equation 

(5)) 

equation (5) as compare! 
with equation (7) 
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TABLE 111.- Od AS A FUNCTION OF MACH NUMBER FOR WAKE FORM IN 

CASE III 


Pl~Po 

^o“Po 

Mach 

number 

. 

c d 

Equation 

(6) 

e a 

Equation 

( 5 ) 

Percentage of error In 
equation ( 5 ) as compared 
with equation (6) 

0.1 

0.4 

0.0603 

0.0621 

3.0 

.1 

• 6 

.0570 

• 05 S 4 

2.5 

.1 

.8 

.0529 

• OB4o 

2.1 

.1 

1.0 

. 04 . 3 B 

. 0490 ' 

1.0 

0 

.4 

.o64g 

.0662 

2.2 

0 

.6 

.0607 

.0617 

1.6 

c 

.S 

.0556 

.0562 

1.1 

0 

1.0 

.0 BOO 

.0601 

0.2 

-.2 

.4 

• 0725 

.0734 

1.2 

-.2 

.6 

.0667 

.0672 

-.7 

-.2 

.S 

• 0593 

.0593 

0 

-.2 

1.0 

.0608 

.. . 0606 .. 

-.6 .. 

-.4 

.4 

.0792 

•0793 

.4 

-.4 

• 6 

.0714 

.0714 

0 

-.4 

.8 

.0612 

.0603 

-.7 

-.4 

1.0 

.0471 

.0464 

- 1-5 


TABLE IV.- cd. AS A FUNCTION OF MACH NUMBER FOR WAKE FORM IN 

CASE IV 


Pi -Po 
^0~Po 

Mach 

number 

c a 

Equation 

(6) 

c a 

Equation 

( 5 ) 

Percentage of error in 
equation (5) as compared 
with equation ( 6 ) 

0 

0.4 

0.0223 

0.0229 

0.4 

0 

. 6 

.0211 

.0212 

• 5 

0 

.8 

.0190 

.0190 

0 

0 

1.0 

.0166 

.0163 

0 

-.1 

*4 

• 0239 

.0240 

.4 

-.1 

, 6 

.0219 

.0220 

.6 

-.1 

.3 

.0195 

.0196 

0 3 

-.1 

1.0 

. 0 l 6 S 

.0168 

0 

-.2 

..4 

.0249 

.0250 

.4 

-.2 

. 6 

.0226 

.0226 

0 

-.2 

.8 

.0193 

.0193 

0 

-.2 

1.0 

.0167 

.0166 

-.6 

- . 4 

.4 

.0267 

.0263 

.4 

-.4 

. 6 

.0233 

.0238 

0 

-.4 

.3 

.0201 

.0200 

— . n 

-.4 

- 

1.0 

.0191 

.ci 50 

-•7 


TABLE V.- cd AS A FUNCTION OF MACH NUMBER FOR WAKE FORM IN 

CASE V 


P1-P0 

Kq-Po 

Mach 

number 

c d 

Eauation 

'(6) 

c d 

Equation 

( 3 ) 

Percentage of error in 
equation (5) as compared 
with equation (6) 

0 

0.4 

0.1046 

0.0967 

" 7-7 

0 

.6 

.0982 
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Figure i - Diagram Of Airfoil, With Wake, In A Compressible Fluid 
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Figure 7. - Types Of Wake Shapes Employed. 
da.) Cases I, II, And EL 




